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Abstract 

This  paper  carries  out  numerical  study  of  the  flow  separations  around  NACA 
0012  airfoil  at  large  angles  of  attack.  Flow  separation  introduces  two  major 
effects:  sudden  loss  of  lift  and  generation  of  aerodynamic  noise.  These  two 
factors  are  highly  concerned  on  the  aircraft  designation.  This  study  gives  a detail 
picture  of  flow  separation. 

As  the  problems  caused  by  flow  separation  are  complicated,  the  spatial  and 
temporal  complexity  makes  it  difficult  to  access  by  conventional  experiment 
methods.  In  the  presented  work,  the  numerical  investigation  results  from  solving 
the  time-dependent  Navier-Stockes  equations  in  the  generalized  curvilinear 
coordinates.  Using  a fourth  order  centered  compact  scheme  for  spatial 
discretization  facilitates  high  resolution  of  the  flow  field,  which  will  be  neglected 
if  using  low-order  numerical  schemes.  To  avoid  possible  non-physical  wave 
reflection,  the  non-reflecting  boundary  conditions  are  used  at  far-field  and  outlet 
boundary. 

Complex  flow  separation,  vortex  shedding,  vortex  merging,  and  vortex  paring 
are  observed  in  the  computational  results.  The  main  purpose  of  this  paper  is  to 
provide  more  detailed  information  of  the  flow  separation. 

1.  Introduction 

Flow  separation  around  airfoil  receives  highly  concerned  from  the  researchers. 
Shih  et  al.(  1 992)  had  pointed  out  that  the  understanding  of  this  problem  needs  to 
make  a step  forward  from  qualitative  conjecture  to  quantitative  measurement  of 
the  instantaneous  flow  field.  Shih  et  al  (1995)  investigated  the  unsteady  flow  past 
a NACA  0012  airfoil  in  pitching-up  motion  experimentally.  They  carried  out  the 
experiment  at  a Reynolds  number  of  5000  using  PIDV  method.  Instantaneous 
velocity  field  data  at  different  times  have  been  obtained  over  the  whole  flow  field. 
They  drew  the  conclusion  that  boundary-layer  separation  near  the  airfoil  leading 
edge  leads  to  the  formation  of  a vortical  structure.  The  evolution  of  the  vortex 
along  the  upper  surface  dominates  the  aerodynamic  performance  of  the  airfoil. 
Complete  stall  emerges  when  the  boundary  layer  near  the  leading  edge  detaches 
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from  the  airfoil,  under  the  influence  of  the  vortex.  Furthermore  in  1995  they 
studied  the  stall  process  using  a water  towing  tank  facility,  and  found  that  near 
the  leading  edge  large  vortical  structures  emerge  as  a consequence  of  Van 
Dommelen  and  Shen  type  separation  and  a local  vorticity  accumulation.  They 
pointed  out  that  the  trailing  edge  only  play  a secondary  role  on  the  dynamic  stall 
process. 

Tenaud  & Phuoc  (1997)  used  large  eddy  simulation  (LES)  to  study  this  problem. 
They  described  that  three  flow  fields  can  be  distinguished  according  to  different 
flow  structures.  Near  the  leading  edge,  vortex  shedding  due  to  separation  of  the 
boundary  layer  is  dominant.  In  the  second  area,  which  is  the  middle  part  of  the 
upper  surface  of  the  airfoil,  the  eddy  structures  grow  and  move  downstream.  The 
last  field  is  close  to  the  trailing  edge,  where  the  alternate  vortices  are  created. 

This  work  focuses  on  numerical  simulation  of  flow  separation  around  a NACA 
0012  airfoil  at  a 16  angle  of  attack.  From  the  previous  understanding  of 
experimental  and  numerical  investigation  of  this  problem,  it  is  known  that  large 
vortices  are  generated  and  shed  from  the  leading  edge,  thus  leads  to  fluctuation 
of  pressure  and  the  lift  lost.  But  the  detail  of  this  process  is  still  unclear,  the  main 
purpose  of  this  work  is  to  give  detailed  picture  of  flow  separation  structures 
using  a high-order  and  high-resolution  approach. 

2.  Basic  Equation  And  Numerical  Methods 

The  governing  equations  are  the  two-dimensional  compressible  Navier-Stocks 
equations  in  the  generalized  curvilinear  £-t|  coordinates,  and  in  a conservative 
form: 

1 3(2  d(E-Ev)  d(F-Fv ) 

J dt  dr\ 
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Where  J is  the  Jacobian  of  the  coordinate  transformation,  and  8,x  i^y,T}x,T]y  are 
coordinate  transformation  metrics,  p is  density,  p is  pressure,  u and  v are 
components  of  velocity.  U = u8)x  +v£  , , V = uT}x  + vljy . e is  the  total  energy. 
The  components  of  viscous  stress  and  neat  flux  are  denoted  by  ,T and 

qx,qv,  respectively. 

The  second  order  Euler  Backward  scheme  is  applied  to  solve  Eq.  (1),  that  is, 

3<2',+l  _4<2n  +Q”1  | d{En+{-Ef)  | d(F"+'  - Fvn+l)  _ Q 
2JAt  3£  drj 

Q"+l  is  estimated  iteratively  as:  Q”+l  =QP  + 8QP  , SQP  —Qp+l  -Qp 
Flux  vectors  are  linearized  by  the  local  Taylor  expansion  about  Q as  following: 

En+l  ~EP  +AP8QP  , Fn+]  ~ Fp  + Bp8 Qp 

Now  equation  (2)  turns  into: 

[-1 + AtJ(D,A  + D„B)]8Qp  =R  (3) 

2 ? n 

Where  R = -Aq " ~ 2Q"  + \ Q"' ) - mD(  (E  - E, ) + D,  (F  - F, )]' 

The  superscript  p stands  for  iteration  step.  and  Di;  are  partial  differential 
operators  in  the  % and  r|  directions,  respectively.  The  right  hand  side  of  Eq.  (3)  is 
discretized  using  the  fourth-order  compact  scheme.  A sixth-order  compact  filter 
is  used  to  depress  the  numerical  oscillation.  The  left  hand  side  of  Eq.  (3)  is 
discretized  following  the  LU-SGS  method  (Yoon  & Kwak,  1992).  Then  the  finite 
difference  expression  of  Eq.  (3)  can  be  written  as: 
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For  subsonic  flow,  u,  v,  T are  prescribed  at  the  upstream  boundary,  p is  obtained 
by  solving  the  modified  N-S  equation  based  on  characteristic  analysis.  On  the 
far-field  and  out-flow  boundary,  the  non-reflecting  boundary  conditions  are 
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applied.  Adiabatic,  non-slipping  condition  is  used  for  the  wall  boundary.  All 
equations  of  boundary  conditions  are  solved  implicitly  with  internal  points. 
Specific  details  of  boundary  treatment  can  be  found  in  Jiang  et  al.  (1999) 

3.  Grid  Generation 

A grid  generation  method  first  proposed  by  Spekreijse  (1995)  is  used  to  generate 
the  C-type  grids.  The  grid  generation  method  is  a composite  method,  which 
consist  of  an  algebraic  transformation  and  an  elliptic  transformation.  The 
algebraic  transformation  maps  the  computational  space  onto  a parameter  space, 
and  the  elliptic  transformation  maps  the  parameter  space  onto  the  physical 
domain  by  solving  a set  of  Poisson  equations.  The  orthogonality  of  grids  on  the 
surface  and  near  the  boundary  is  achieved  by  re-configuration  of  the  algebraic 
transformation.  The  grid  numbers  are  841  in  the  £>  direction,  and  141  in  the  T] 
direction.  The  overview  of  the  C-grid  and  grid  near  the  airfoil  surface  are 
displayed  in  Figure  1(a)  and  (b). 

4.  Computational  Result 

The  flow  field  around  an  NACA-0012  airfoil  is  analyzed  by  solving  the  Navier- 
Stocks  equations  using  the  finite  difference  scheme.  Flow  separation  at  a large 
angle  of  attack  (a  =16°)  has  been  studied  using  a high-resolution  numerical 
simulation.  In  this  case,  the  fluid  flow  around  the  airfoil  becomes  very  unstable 
and  different  eddy  structures  are  formed  in  the  vicinity  of  the  airfoil.  The 
Reynolds  number  based  on  the  chord  length  and  the  free-stream  velocity  is  5x10  s, 
the  free  stream  Mach  number  is  0.4.  The  angle  of  attack  is  16  °.  During  the 
computation,  the  flow  field  is  recorded  every  1000  time  steps,  and  the  time  step 
is  approximately  1.768 xlO^L/C/^  . Figure  2 shows  the  contours  of  the 
instantaneous  span  wise  voriticity.  The  time  interval  between  each  of  those  15 
pictures  is  about  0.1012L/U „ . From  those  figures  we  can  see  that  flow 
separation  process  start  at  the  leading  edge  and  the  leading  edge  vortices 
continue  to  shed  and  convect  downstream.  In  Figure  2(c),  a separation  bubble 
can  be  observed  on  the  upper  surface  of  the  airfoil  near  the  leading  edge.  A chain 
of  vortical  structures  appear  on  the  upper  surface  of  the  airfoil,  which  is  more 
clear  in  Figure  2(h).  The  vortices  shedding  from  the  leading-edge  rotate  in  the 
clockwise  direction.  Near  the  surface  of  the  airfoil,  a layer  of  reversed  vorticity  is 
induced  by  the  vortices  shedding  from  the  leading-edge.  The  interactions 
between  the  positive  and  negative  vorticity  leads  to  the  vortex  pairing.  The 
merging  of  vortices  rotating  in  the  same  direction  is  also  observed,  e.g.  in  Figure 
2(j).  A vortex,  which  first  appears  at  the  trailing  edge  of  airfoil  in  Figure  2(a), 
grows  as  it  is  carried  downstream  by  the  mean  flow,  as  shown  in  Figure  2(b),  (c), 
and  (d).  In  Figure  2(d),  along  the  shear  layer  starting  from  the  trailing  edge,  a 
series  of  small  vortical  structures  are  generated  as  a result  of  the  Kelvin- 
Helmholtz  instability.  These  small  structures  are  also  visible  in  Figure  2(e),  (f), 
(g),  and  (h). 
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5.  Conclusions 

Detailed  numerical  simulation  has  been  carried  out  by  solving  Navier-Stocks 
equations  in  the  generalized  curvilinear  coordinates  to  study  the  separated  flow 
around  an  NACA  0012  airfoil  at  large  angle  of  attack.  By  using  a fourth-order 
centered  compact  scheme  for  spatial  discretization,  the  small-scale  vortical 
structures  are  resolved,  which  will  dissipate  if  low-order  numerical  schemes  are 
used.  Non-reflecting  boundary  conditions  are  imposed  at  the  far  field  and  outlet 
boundaries  to  avoid  possible  non-physical  wave  reflection. 

The  numerical  simulation  results  clearly  describe  the  flow  separation  process  at 
the  upper  surface  of  the  airfoil.  The  phenomena  of  the  leading-edge  separation, 
vortex  shedding,  vortex  merging,  vortex  pairing,  and  formation  and  shedding  of 
large-scale  trailing  edge  vortex  are  displayed  and  discussed  in  detail.  The  small- 
scale  vortices  associated  with  the  Kelvin-Helmholtz  instability  are  also  observed 
along  the  shear  layer  near  the  trailing-edge.  These  phenomena  are  in  food 
agreement  with  the  experimental  results  obtained  by  Shih,  et  al  (1992,1995)  and 
Yoshifumi  et  al  (1986). 
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(a)  overview  of  the  grids 


(b)  grids  near  airfoil  surface 


Figure  1.  C-grid  around  an  NACA  0012  airfoil 
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Figure  2.  Contours  of  instantaneous  spanwise  vorticity  at  different  time 


